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We consider the finite-size corrections in the Dicke model and determine the scaling exponents 
at the critical point for several quantities such as the ground state energy or the gap. Therefore, 
we use the Holstein-PrimakofT representation of the angular momentum and introduce a nonlinear 
transformation to diagonalize the Hamiltonian in the normal phase. As already observed in several 
systems, these corrections turn out to be singular at the transition point and thus lead to nontrivial 
exponents. We show that for the atomic observables, these exponents are the same as in the Lipkin- 
Meshkov-Glick model, in agreement with numerical results. We also investigate the behavior of 
the order parameter related to the radiation mode and show that it is driven by the same scaling 
variable as the atomic one. 
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Superradiance is the collective decay of an excited pop- 
ulation of atoms via spontaneous emission of photons. 
This phenomenon first predicted by Dicke in 1954 Q has, 
since then, been observed experimentally in several quan- 
tum optical as well as solid-state systems (for a review 
see Ref. 0)- The phase diagram of the Dicke model, 
which is the subject of the present study, has been estab- 
lished in the thermodynamical limit by Hepp and Lieb Q 
revealing the existence of a second-order quantum phase 
transition. This transition has been shown to be associ- 
ated to a crossover between Poisson and Wigner-Dyson 
level statistics for a finite number of atoms N , thus rais- 
ing the question of the finite-size corrections in this sys- 
tem 0, l^l ■ These corrections have also been shown to be 
crucial in the understanding of entanglement properties 

01 which become trivial if one directly considers the 
thermodynamical limit 0,0- In these latter studies, non- 
trivial finite-size scaling exponents have been numerically 
found at the critical point and further been compared to 
those obtained in the Lipkin-Meshkov-Glick model 
The aim of the present work is to determine these expo- 
nents. 

To achieve this goal, we proceed in several steps. First, 
we use the Holstein-Primakoff boson representation 
for the atomic degrees of freedom which is well adapted 
for a 1/N expansion of the Hamiltonian, N being the 
number of atoms. Second, we exactly diagonalize the 
expanded (quartic) Hamiltonian at order In a re- 

cent series of papers 0, 0, this diagonalization 
was performed using the Continuous Unitary Transfor- 
mations (CUTs) methods but here, the problem is 
more complicated for several reasons: (z) it involves two 
different degrees of freedom; {it) the parameter space is 
two-dimensional; and (Hi) the total number of particle 
is not fixed. These complications render the analytical 
resolution of the flow equations coming from CUTs ap- 
proach difficult ^3 ■ We are thus led to use an alternative 
approach relying on a canonical transformation of the 



initial bosonic operators. This transformation provides 
both the eigenstates and the eigenspectrum of H, and 
thus allows one to compute any matrix element of any 
observable. Here, we focus on the quantities which have 
been numerically investigated and we show that their 
expansion is singular at the critical point. The anal- 
ysis of these divergences directly provides the finite-size 
scaling exponents which are the same as in the Lipkin- 
Meshkov-Glick model, at least for the physical quantities 
involving atomic degrees of freedom. We also compute 
this exponent for the order parameter which is found to 
vanish as N~^^^ at the transition point. Finally, we dis- 
cuss numerical data which are in good agreement with 
our predictions. 

Let us consider the single-mode Dicke Hamiltonian Q 
without the rotating wave approximation 

H = luqJ^ + Lua'^a + ^= {a'' + a) {J+ + J^) , (1) 



/2j 

where a' and a are bosonic creation and annihilation 
operators satifying [a,a^] = 1. The angular momentum 
operators are defined as Ja = J2iLi '^a/'^ where the CTa's 



are the Pauli matrices, and J± 



This Hamiltonian, which describes the interaction of 
a photon field with N two-level atoms (spins 1/2), con- 
serves the magnitude j of the pseudo-spin ( [H, J^] = O) . 
In the following, we focus on the sector j ~ N/2 to which 
the ground state belongs. Further, one has [-ff, H] = 
where 



H : 



(2) 



is the parity operator. An appropriate basis of the 
Hilbert space is thus provided by the states \n) (g) \ j,m) 
where \n) denotes an eigenstate of the photon density op- 
erator a^a with eigenvalue n, and \ j, m) the eigenstate of 
and Jz associated to eigenvalues j and m respectively. 
In the thermodynamical limit and at zero temperature, 
the system described by this Hamiltonian undergoes a 
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second-order quantum phase transition at a critical cou- 
plings Ac = As an order parameter of the tran- 
sition, one can choose the expectation value of the photon 
number per atom in the ground state which satisfies: 



independent coefficients. Indeed, one has schematically: 



N 



lim {a) a) IN = 



for A < Ac 
_ _^ for A > Ac 



(3) 



As we shall see, nontrivial exponents are only found at 
the critical point that we will investigate from the normal 
(symmetric) phase, i.e., for A < Ac- A convenient start- 
ing point to perform a \/N expansion of the Hamiltonian 
is to use the Holstein-Primakoff boson representation of 
the angular momentum which reads: 



J+ = bWN -b^b = {J^)\ 
N 

Y' 



J. 



6^6 



(4) 
(5) 



with [b, b^] = 1, so that we now have to consider a two- 
boson problem. In the thermodynamical limit and for 
A < Ac, one has {b^b)/N ^ 1 and we can expand the 
square root in Q to obtain the following expanded form 
of the Hamiltonian: 



H 



N 

'2N 



ujo + LUob^'b + oja^a + X (a^ + a) (5^ + b) 
{a^ +a) (^b^b^ + b^%'^ +0{l/N^) . (6) 



Note that we restrict this expansion at the order 
which, as we will see thereafter, is sufficient for our pur- 
pose. At order (1/A^)°, the Hamiltonian is quadratic and 
can thus be diagonalized via a Bogoliubov transforma- 
tion as already discussed in Ref. [5|. The real problem 
arises at the order where one has to diagonalize a 
quartic form. 

As explained above, the CUTs formalism used in recent 
studies [l^ |lj for this step is difficult to implement 
in the Dicke model. Instead, we use here an approach 
that simply requires to solve a set of algebraic equations 
instead of differential equations. The main idea of this 
method is to perform the following canonical transforma- 
tion 



A 



6t 



T — 



(7) 
(8) 



where the A^^ and i?^ are polynomials functions of new 
bosonic operators c^,c,c?^,d, such that H expanded at 
order l/N^ is a polynomial function in Uc and rid- 

At order zero, this transformation coincides with the 
Bogoliubov transformation and one has to determine 8 



4 = E -^L^'^^'d'^ 

i,j,k,l 
i,j,k,l 



(9) 
(10) 



where ctf^ji^i (resp. Pl'^J i) stands for the coefficient of 

c^^c'd'^^d} in the expansion of A]^ (resp. i?J). Since, 
at this order, the transformation is linear, the sum is 
constrained by z + j + fc + Z = 1. The eight equations to 
be solved which are quadratic forms of the a' s and /3's 
are, as usual, obtained by (i) requiring the cancellation 
of (nonconstant) terms which are not proportionnal to 
TT-c and rid, and (ii) imposing the following commutation 
rules, 

[a, a)] = 1, [b, b^] = 1, [a, b^] = 0, [a, 6] = 0. (11) 

The full solution of these equations can be found in 
Ref. 0. 

Now, let us turn to the next order p = 1 for which H is 
quartic. At this order, the corresponding transformation 
reads 



4 = E "Sm^^'^^^'^'' 

i,j.k.l 



(12) 
(13) 



where the sum now contains two types of terms: linear 
{i+i+k+l — 1) and cubic {i+j+k+l = 3). There is thus 
48 independent parameters to be determined. At this 
order, these are the only terms that need to be present 
since the Hamiltonian only contains quadratic and 
quartic terms. We also emphasize that once the a| ;'s 
and the (3^^ ;'s are known, the constraints to be satisfied 

arc linear functions of the ct-^^-^^iS and /?! ^''j, ;'s. More 
generally, to determine the parameters for p> 1 , we must 



solve a set of linear equations involving only the af^^ ^ ; 's 
with q < p. At order p = 1, the equations to be solved 
are given by requiring the cancellation of (nonconstant) 
terms not proportionnal to ric, Ud, n^, n^, and ncUd, 
but also by requiring the commutation rules (|ll|l to be 
satisfied. Note that the spirit of this approach is the same 
as the one issued from the CUTs in which the running 
coupling, in the infinite time limit, identify with the a 
and /3's lli. 

The exact solutions of this set of equations are ob- 
viously too long to be given here, but let us sketch the 
main results that can be extracted from them. As already 
shown in several models 

113,111113: the 1/A^ corrections 
to physical observables such as the gap or the order pa- 
rameter display some singularities at the critical point. 
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As detailed in 0|, the schematic form of an observable 
$ in the vicinity of the critical point is: 



$A.(A) = <i>"^^(A) + <i>?7s 



(A), 



(14) 



where the superscript reg and sing stands for regular and 
singular functions at A = Ac- By singular, we mean that 
the function and/or its derivatives with respect to A di- 
verges at the critical point. Further, a close inspection of 
the expansion shows that near Ac one has: 



$7^(A) 



S(A)«* 



iVS(A)3/2 



(15) 



where E!(A) = Ac — A and is a function depending only 
on the scaling variable N'E{X)^/^. The exponents ^$ and 
n$ are characteristics of the observablcs $. In the present 
study, we have only checked this scaling hypothesis at 
order 1/A^ but we strongly believe that, as in previous 
models we studied, one indeed has such a scaling variable. 
For instance, the ground state energy per atom near the 
critical point reads: 



eo 



ci + C2S(A)i/2 



1 C3 

iV2 S(A) 



(16) 



with 



Co 

Cl 
C2 

C3 



~^o/2, 



(C^2+ ^2)1/2' 

3a;5/2^g/2 
64(^2 +a;2)' 



(17) 
(18) 

(19) 
(20) 



Using the hypothesis (|15|l . these expressions allow us to 
identify = 1/2 and = 1. Note that for the spec- 
trum (only), one can also obtain these corrections by a 
standard first-order perturbation theory. The most strik- 
ing result is that the scaling variable A^S(A)^''2^ which 
is the key ingredient for our study, does not depend on 
the observable. This remarkable fact already observed 
for single-boson model |lo| is rather surprising here since 
one may have expected one different variable for each 
types of degrees of freedom. Furthermore, the Hamilto- 
nian depends on two independent parameters but their 
value do not change the scaling variable. In particular, 
we find no difference between the resonant {lo = cjq) a-nd 
the off-resonant case. 

To obtain the finite-size scaling exponent from the 
general form (|15|l it is sufficient to underline that, at 
finite N , no divergence can occur in the behavior of the 
observables, even at the critical point. This straightfor- 
wardly implies that, to cure the singularity coming from 
S(A)^*, one must have J^$(x) ^ x'^C^/s^ This behav- 
ior of T then leads to ^'^"'^(xc) - jv-("*+2«*/3). We 



have computed the finite-size scaling exponents for sev- 
eral observables which are summarized in Table ^ For 
completeness, we also give the value of these quantities 
in the thermodynamical limit. 
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TABLE I: Finite-size scaling exponents at the critical point 
for the ground state energy eo, the gap A, the order parameter 
{a^a)/N. the magnetization per atom {J^)/N, and the two- 
point correlation function {Ja)/N'^ for a = x,y,z. 



It is clear that the canonical transformations (|12I13|) 
we used to diagonalize the Hamiltonian at order 1/A^ 
allow us to compute any matrix element (not only diag- 
onal) of any observable expressed in terms of the initial 
operators. Here, we only focused on ground state ex- 
pectation values (except for the gap) because these have 
already been numerically computed and can thus be di- 
rectly checked. 

The finite-size scaling exponents at the critical point 
have been computed for three quantities Q: {Jz)/N 
(-0.54 ± 0.01), y/{Si)/N (-0.35 ± 0.01) and indirectly 
{Jy)/N (—0.26 ±0.01). These results are very close to our 
predictions which are —2/3, —1/3 and —1/3 respectively, 
as can be read in Tabled Nevertheless, it is true that 
our results do not lie within the error bars proposed by 
Reslen et ai. The same discrepancy was already observed 
in the LMG model for which we have explicitely shown 
that it was due to the too small system sizes investigated 
[l(l H^. Here, we strongly believe that the asymptotic 
regime was also not reached but, unfortunately, it is diffi- 
cult to consider significantly larger sizes as those studied 
in Ref. |^ . This clearly requires further numerical efforts 
[l5l | which are beyond the scope of the present study. 

Let us also mention that the concurrence C studied in 
Ref. which measures the spin-spin entanglement [l^ 
reads 



{N -1)C =1-4{J^)/N. 



(21) 



We thus predict a finite-size scaling exponent for this 
(rescaled) concurrence which is —1/3. 

At first glance, these results are strikingly similar to 
those obtained in the LMG model 0, Q| and this calls 
for several comments. Indeed, it is well-known that if one 
focuses on the atomic degrees of freedom, both systems 
are equivalent in the thermodynamical limit as shown 
with different methods |^, 0, lisl Il9| . However, the 
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finite-size corrections fail to be captured through this 
mapping. For instance, in the Dicke model, one has 
limTv^oo 4( ) /N = ujnl(ijp- + i^q)"^^^ whereas it vanishes 
in the LMG model [la, Q. Moreover, for the LMG 
model, these exponents were found to be related to the 
upper critical dimension and the mean-field critical expo- 
nents of the Ising model in a transverse magnetic field |20j 
which is the counterpart of the LMG model with short- 
range interactions. For the Dicke model, it is difficult to 
find such a mapping since one cannot simply consider it 
as a long-range interacting system which would admit a 
short-range equivalent. Consequently, the similarity be- 
tween the exponents of these two models is a nontrivial 
result which shed light on a recent controverse on that 
subject 

Unlike previous studies using CUTs, we have devel- 
oped here an alternative simple perturbative approach 
relying on a canonical transformation which allows one to 
diagonalize the Hamiltonian at order 1/iV. This method 
can, in principle, be applied to many similar models in- 
volving more than one type of boson and requires to solve 
a set of linear equations. It is thus, a 'priori simpler than 
the CUTs technique even if the number of equations to 
be solved quickly grows with the order of the ex- 
pansion. Whatever the approach chosen, the main result 
to keep in mind is that if one accepts the hypothesis of 
a unique scaling variable, it is sufficient to compute the 
first nontrivial correction of one observable (for example 
the ground state energy) to get all the exponents. In- 
deed, the determination of and for the other ones 
can already be infered from the quadratic approximation. 

Finally, let us quote a recent work |24j where a semi- 
classical approach has been introduced to obtain the 
finite-size scaling exponent in the LMG model. It would 
be interesting to analyze the Dicke model within this 
framework to have a better understanding of the simi- 
larities between these two systems. 
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